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A constitutive theory for use in structural and durability analyses of 
high-temperature isotropic alloys is presented. Constitutive equations based 
upon a potential function are determined from conditions of stability and phys- 
ical considerations. The theory is self-consistent; terms are not added in an 
ad hoc manner. It extends a proven viscoplastic model by introducing the 
Kachanov-Rabotnov concept of net stress. Material degradation and inelastic 
deformation are unified; they evolve simultaneously and interactively. Both 
isotropic hardening and material degradation evolve with dissipated work which 
is the sum of inelastic work and internal work. Internal work is a continuum 
measure of the stored free energy resulting from inelastic deformation. 


INTRODUCTION 

The nucleation, growth and coalescence of voids and microcracks are physi- 
cal phenomena that degrade a material's continuity. This degradation results 
in a loss of strength, and is the eventual cause of failure. Continuous damage 
mechanics applies whenever the distribution of defects does not include one or 
more dominating macroscopic cracks; otherwise, fracture mechanics applies. The 
subject of this paper falls under the topic of continuous damage mechanics; 
applications to fracture mechanics are not discussed. 

A constitutive theory applicable to structural and durability analyses of 
high-temperature Isotropic alloys is developed. A set of constitutive equa- 
tions based on a single potential function is determined from stability con- 
ditions and physical considerations. A specific potential function from a 
proven viscoplastic theory is extended to account for internal damage by intro- 
ducing the Kachanov - Rabotnov (refs. 1 and 2) concept of a net stress. Inter- 
nal damage and inelastic deformation are unified in this approach; they evolve 
simultaneously and interactively. The theory is self-consistent in that it is 
derived from a potential function; terms are not added in an ad hoc manner. 
Other viscoplastic theories that incorporate continuous damage mechanics have 
been proposed. The evolutionary equations for material degradation in the 
theories of Chaboche (ref. 3), Bodner (ref. 4), and Walker and Wilson (ref. 5) 
are phenomonological ly determined, whereas, the Perzyna theory (ref. 6) is 
micromechanistical 1y based. In this paper the evolutionary equation for mate- 
rial degradation is derived from a potential function. 


‘NASA Lewis Resident Research Associate. 
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Continuous damage mechanics deals with, an infinitesimal volume element 
of material (called a particle), whose dimensions are large enough to contain 
many material defects, yet small enough to be treated as a mathematical point 
in a continuum sense. Consider a face on such an element whose unit normal is 
given by n^.^ Let A denote its surface area in a flawless (or undamaged) 
state, and let A' denote its net surface area in the presence of material 
defects (or in a damaged state); thus A' < A. The internal damage associated 
with this particle, and in the direction of this unit normal, is defined by 


(i) 


A - A’ 


( 1 ) 


which is bounded by the interval 0 < w < 1 where u » 0 in an undamaged 
state. Whenever the orientations of material defects have preferred direc- 
tions, damage becomes a function of these directions resulting in an entity of 
tensorial nature (refs. 7 to 9); otherwise, damage is isotropic and can be 
represented by a scalar. In this paper, damage is taken to be isotropic as a 
simplifying assumption. Kachanov (ref. 1) calls the quantity «|» « 1 - u the 
continuity of the material. 


Consider once again a face on an infinitesimal material volume element. 

In an undamaged state, traction is the ratio of the force transmitted through 
the surface to the surface area A. It is related to the unit normal n; 
by a homogeneous linear operator a,, called the applied (or Cauchy) stress, 
that is2 


i 




( 2 ) 


In a damaged state, traction becomes the ratio of the force transmitted through 
the surface to the net surface area A'. It is related to the unit normal 
n^ by a homogeneous linear operator called the net (or Kachanov- 

Rabotnov) stress; thus 


F 


1 


A' 




Combining equations (1) to (3) results in 


(3) 


ij * (1 - 


which relates the net stress to the applied stress. 

Like the classical theories of creep and plasticity, strain 
given by the sum 'J 


(4) 


Ull scalar, vector and tensor fields are defined at particles whose spa- 
tial coordinates are at the instant t in a Cartesian reference frame. 

^Repeated indices are summed over in the usual manner. 
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(5) 


e p 

where is the elastic strain and is the inelastic (or plastic) 

strain. From a thermodynamic viewpoint, the elastic change in strain is the 
reversible portion of a change in strain, while the inelastic change in strain 
is the irreversible portion of that change in strain. Small displacements and 
rotations are assumed. 


ship 


The elastic response of polycrystalline metals is given by the relation- 


■ij 


<1 . . « AT 

Ed - (o) 


( 6 ) 


where v is the Poisson ration, a is the mean coefficient of thermal expan- 
sion, AT is the temperature change and is Kronecker's delta. Since E 

is Young's modulus in an undamaged state, and E' = Ed - «) can be considered 
as Young's modulus in a damaged state, we obtain the following expression: 


(E - E') 
“ E 


(7) 


This is a useful measure of internal damage, because it can be readily deter- 
mined by experiment (ref. 10). 


CONSTITUTIVE THEORY 

Much of the essential structure in the classical theory of plasticity 
derives, not so much from thermodynamic concepts, but from concepts of mate- 
rial stability as described by Drucker (ref. 11). A single postulate of 
stability is sufficient to unify the description of inelastic behavior of time- 
independent materials under isothermal conditions. A dual postulate of stabil- 
ity has been applied by Ponter (ref. 12) to time-dependent materials whose 
hereditary behavior can be represented in terms of internal state variables Ea 
(a a l,2,...,n) and their conjugate thermodynamic forces f^. In that work, 
small isothermal changes in stress at constant internal state are assumed to 
obey the inequality 


do.^ dej. > 0 (8) 

where f and T are constant; whereas small isothermal changes in internal 
state at constant stress are assumed to satisfy the inequality 

df dE >0 (9) 

a a - 

where and T are constant. In contrast, a thermodynamic counterpart 

to the second-order inequality in equation (9) is the restriction of positive 
internal dissipation 
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(10) 


f ^ >0 

a^a - 

which Is derived from the second law. 

Since changes in inelastic strain rate and internal variable rates are 
path independent in the complete state space, the inequalities in equations (8) 
and (9) can be integrated along a straight-line path between two arbitrary 
11 2 2 

states (a^ .f ,T) and <a_,f ,T) resulting in the following inequality: 
ij a ij a 




(11) 


Along a constant stress path (i.e. under conditions of creep) only the last 
term in this inequality remains, and we can easily show that a sufficient con- 
dition for its satisfaction is 


1 90 

“ af 


( 12 ) 


a 


where 0(a|^ ,f< 3 ^,T) is convex and positive definite in a.j,f< 3 j. Here we assumed 

that equation (12) is not constrained just to constant stress conditions, but 
is valid in general. Rice (ref. 13), Martin (ref. 14), and others have shown, 
using thermodynamic arguments, that if the kinetic (or evolutionary) law can be 
expressed as equation (12), then the flow law given by 


-P _ 9Q - 
'U - 


(13) 


is a derived result. 


The criteria for stability and the resulting kinetic and flow laws lead to 
a vital theorem (ref. 12): 


"The stress and state histories are uniquely defined for time t > to by 
the initial conditions at t = t© and the loading history." 

The existence of this theorem is essential if this is to be a meaningful con- 
stitutive theory for use in structural analyses. 

Following the lead of Ponter and Leckie (ref. 15) and Ponter (ref. 12), we 
adopted an additional constitutive assumption, that is 


g : ^ 

h(f ) “ “^g ' " 8f 

g g 


(14) 


in which h is a hardening function of the Internal force f^. The physical 
origin of equation (14) in describing the local response of a crystallographic 
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slip system, and the limitations that result In transferring from a local 
format to a global one, are discussed by Ponter and Leckle. An additional 
motive for adopting equation (14) comes from considering conditions In the 
neighborhood of a stress free state, as depicted In figure 1. In particular, 
for a "J 2 -type" material (considered In the following section), the surface 
(0(0,fa,T) ■ constant) Is a sphere of radius If^l In thermodynamic force 
space. The gradient vector 80/8fa at each point on the surface where 0 - 
constant Is directed along the outward normal. By considering the constitu- 
tive assumption fgj/h - the thermodynamic restriction In equation (10) 
can be expressed as fa^a^t^ 1 0, which for positive h constrains the vector 

to be contained within a half-sphere In thermodynamic force space. (See 
ng. 1.) The Ponter-Leckle constitutive assumption equation (14) selects the 

direction of f<x/h so that Its projection on f<x Is a maximum; that Is, the 
Ponter-Leckle constitutive assumption ensures that state recovery occurs under 
maximum Internal dissipation In the neighborhood of a stress free state. 

The extended normality structure expressed In equations (13) and (14) 
provides the basis for the present development. Moreover, this structure is 
assumed to hold under noni sothermal conditions. 


A SPECIAL POTENTIAL FUNCTION 

The governing differential equations of a theory of viscoplasticity that 
accounts for Internal damage are taken to be associated with the normality 
structure of a potential function Q as discussed in the previous section. 

The independent arguments of this potential function are the applied stress 
oM, an internal stress (3)j, a threshold strength Z, the damage u, and the 
temperature T; thus, n(ai j,(ii j ,Z,u) , where the temperature dependence 1s 
Implicit. From a thermodynamic viewpoint, the Internal stress and the thresh- 
old strength are averaged thermodynamic forces, and damage Is an averaged 
Internal variable (or thermodynamic displacement). The Internal stress and 
threshold strength are associated with kinematic and Isotropic hardening 
behaviors, whereas damage Is associated with material degradation. 

Moderate states of hydrostatic pressure have virtually no Influence on 
the Inelastic response of metals. ^ The stress dependence of Q can therefore 
be expressed in terms of the deviatoric applied stress 

^ij “ ''ij ■ I "'kkSj 
and the deviatoric Internal stress 

®ij “ l^lj ■ i l^kk^lj 


(15) 


(16) 


^Moderate states of hydrostatic pressure have a strong Influence on the 
formation and growth of material defects and, therefore, on the damage. This 
effect, however. Is accounted for in the degradation function, not in the 
potential function. 
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where their difference 


is the effective stress associated with inelastic deformation. Under condi- 
tions of full isotropy, the invariants 

^2 “ i ®ij®ji 

(18) 

^3 “ I ®l3^jk®ki 
and 

^2 " 1 ^ij^ji 

(19) 


'^3 ■ 3 ^ij^jk^ki 

provide a complete description of the stress dependence of Q. 

For the chosen potential function, the flow law is 

( 20 ) 


( 21 ) 
(22) 

(23) 

in accordance with the results of the previous section. Here h^ and h 2 
are the hardening functions for the internal stress and the threshold strength, 
and D is a degradation function. Equation (20) is the flow law of Rice 
(ref. 13). Equations (21) and (22) are the evolutionary laws of Ponter and 
Leckie (ref. 15) and equation (23) is the proposed evolutionary law for damage. 


'13 - 

and the evolutionary laws are taken to be 


^ii an 

L ^ 

h^ * ■ az 
and 

D ” aw 
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since damage is an internal variable, whereas Internal stress and thresh- 
old strength are thermodynamic forces, the sign for equation (23) is different 
from that for equations (21) and (22). The reason for this difference is a 
Legendre transformation, like those used in equilibrium thermodynamics. These 
equations form the foundation for a theory of viscoplasticity that Incorporates 
internal damage. A specific model is obtained by choosing a particular form 
for the potential function. 

The potential function considered for this model is 

£2 - Jk^ f(F) dF + Jk^ g(G) dG + Jz(Z) dZ (24) 

where the stress dependence enters through the functions F(Zij) and G(Bij). 
This extends the function used by Robinson (ref. 16) to include Isotropic 
effects. The fact that equation (24) is a sum of integrals is consistent with 
Rice's formulation (ref. 13). In his definition of the potential function, 
each integrand denotes the rate of change of a thermodynamic displacement (or 
internal variable), which is Integrated with respect to its conjugate thermo- 
dynamic force. 

In the spirit of von Mises (ref. 17), the stress dependence of F and G 
relies only on the second Invariants; in particular,^ 

F - -f - 1 (25) 

and 





Equation (25) is a Bi ngham-Prager (refs. 18 and 19) yield condition with K 
denoting the yield strength in shear. Inelastic strain only occurs when 
F > 0; an elastic domain is defined by the inequality F < 0. The boundary 
between these two regions, F ■ 0, is a sphere in deviatoric stress space; it 
is the threshold or quasi-static yield surface. The origin of this sphere is 
at Bm, and its radius is K. The inelastic domain, at a fixed inelastic 
state, '^consists of a nested family of spherical surfaces in deviatoric stress 
space; each is a surface of constant F, and thus of constant £2. Viscoplas- 
ticity differs from classical plasticity in that stress states that lie out- 
side the quasi-static yield surface are admissible; they are not admissible in 
classical plasticity. 


^Many theories of viscoplasticity take F « J2/<2 Instead of equa- 
tion (25); thus, there Is no elastic domain. The only influence that this 
choice for F would have on the resulting theory is that the inequality would 
be removed from the flow function. 
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If the stress dependence of F and G Is to be a net stress dependence 
in the sense of Kachanov (ref. 1) and Rabotnov (ref. 2), then K must be a 
linear homogeneous function of damage, 

K = k(Z)(1 - «) (27) 

where k is the quasi-static yield strength in shear. Although damage, 
by definition, influences the state of stress, this influence is manifested 
by a reduction in strength as exemplified in equations (7) and (27). In 
equation (27) the reduction in strength due to material degradation competes 
with the process of hardening which enhances strength. 

Given the potential function (eq. (24)), the flow law of equation (20) 
becomes Prager's flow equation (ref. 20), 

2p cPj » f(F)Z^^ (28) 

where p is the viscosity and f is the flow function. (This is derived in 
the appendix.) The Bingham-Prager yield condition (eq. (25)), constrains the 
flow function so that it is zero in the elastic domain. Coaxiality between 
the effective stress and the inelastic strain rate is implied in equation (28). 
Stability (in the sense of eq. (8)) constrains the flow function f(F) to be 
nondecreasing with increasing values of F. Most theories of viscoplasticity 
use the general form of this flow equation. 

Given the potential function (eq. (24)) the evolutionary law for internal 
stress (eq. (21)), becomes a Bai ley-Orowan type relationship (refs. 21 and 22), 
that is 


B^^ = 2hjj(G)cfj - r^(G)B^j (29) 

where hj) and r^ are the kinematic functions for hardening and thermal 
recovery. (Equation (29) is derived in the appendix.) The first term in this 
equation, for constant h^, is Prager's rule for kinematic hardening (ref. 20). 
To model dynamic recovery of the internal stress, Robinson (refs. 23 and 24) 
presents a kinematic hardening function that exhibits an analytical discontinu- 
ity whenever there is a reversal in stress. 

The second term in equation (29) accounts for the thermal recovery of the 
internal stress state. This is an anelastic response since it continues until 
the internal stress has relaxed to zero, regardless of whether the current 
deformation state is elastic or inelastic. 

The experimental results of Mitra and McLean (ref. 25) verify the 
Bai ley-Orowan hypothesis that inelastic deformation occurs as a result of two 
competing mechanisms: a hardening process that progresses with inelastic 

deformation, and a thermal recovery process that progresses with time. When- 
ever these two mechanisms balance such that Bij » 0, the internal stress is 
in a steady state. Stability (in the sense of eq. (9)) constrains the func- 
tion g(G) - r|j(G)/2h5(G) to be nondecreasing with Increasing values of G. 
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The general form of this evolutionary equation for the internal stress is used 
in many viscoplastic models. 

Data from metals that strain-age indicate that the evolution of the quasi- 
static yield strength (defined in eq. (27)) depends on the history of thermo- 
mechanical loading (ref. 26). These data suggest an evolution such that 

K = rZ - 0(Z)t (30) 


where the parameter r reflects the temperature dependence of the quasi-static 
yield strength in an annealed state, and the function 0 represents the change 
in quasi-static yield strength resulting from a change in temperature. The 
evolution of the quasi-static yield strength given by equation (30) is path 
independent whenever the following equation is satisfied: 

0(Z) = -Z (31) 


If equation (30) is to be self-consistent, then k must be path independent 
in the annealed state; therefore 


0(Z ) 

d 



(32) 


where 0 < Z^ < Z. This constraint must always be satisfied; it is like an 
initial condition for the functional dependence of 0. 

Given the potential function (eq. (24)) and the equation of evolution for 
the quasi-static yield strength (eq. (30)), the evolutionary law for the thres- 
hold strength (eq. (22)) becomes a Bai ley-Orowan type relationship (refs. 21 
and 22), 


where 


Z = 


h^(Z) 





*^i.i^.ii 
“ 2hb 


(G) 


(33) 


(34) 


and h 2 and r 2 are the isotropic functions for hardening and thermal 
recovery. (These equations are derived in the appendix.) Equation (33) 
implies that the path of thermomechanical loading influences the rate of iso- 
tropic hardening; under isothermal conditions it reduces to 


. h,(Z)W 
Z = -^-r^(Z) 


(35) 
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The first term in equations (33) and (35) implies that isotropic hardening 
progresses with dissipated work. This dissipated work, as defined in equa- 
tion (34), is the sum of the inelastic work and internal work. 5 The internal 
work can be thought of as a continuum measure of the free energy stored in the 
material that arises from inelastic deformation. With the exception of the 
viscoplastic model of Bodner and Partom (ref. 27) (where kinematic hardening is 
not present, and isotropic hardening evolves with inelastic work), all visco- 
plastic models that incorporate isotropic hardening, to the best of our knowl- 
edge, assume that this process progresses with inelastic path length 

J(c^jc^i)^’ dt. This is an assumption that our theoretical derivation does not 
support. 

The second term in equations (33) and (35) accounts for the anelastic 
thermal recovery (or annealing) of the threshold strength. This function must 
be constrained so that recovery terminates when the annealed value of thresh- 
old strength is obtained. Stability (in the sense of eq. (9)) constrains the 
function z(Z) » r 2 (Z)/h 2 (Z) to be nondecreasing with increasing valuei of Z. 

Given the potential function (eq. (24)), the evolutionary law for damage 
(eq. (23)) becomes 


0) - D(w)i-!^*— (36) 

I - (0 

where D is the degradation function. (Equation (36) is derived in the 
appendix.) 

Since materials do not degrade in states of sufficient hydrostatic com- 
pression, in general the degradation function ought to switch off the evolu- 
tion of internal damage when a critical state of hydrostatic compression is 
reached. Equation (36) implies that damage evolves with dissipated work, as 
defined in equation (34). Hereditary effects are included through the depend- 
ence of dissipated work on inelastic strain and internal stress: thus, equa- 
tion (36) has the potential to account for time-dependent effects in a natural 
way. Stability (in the sense of eq. (9)) is satisfied if D(w) does not 
increase with increasing values of «; but this is not observed. Initially 
the dissipation function is virtually a constant, and the material response is 
stable for all practical purposes. However, near the end of life, the value 
of the dissipation function explodes, thereby leading to material instability 
or failure. This is not to say that this theory is undesirable, for it is 
precisely this instability that continuum damage mechanics attempts to 
characterize. 

Many researchers have used inelastic work as a parameter to characterize 
fatigue damage (e.g. refs. 28 to 31). The equation of damage evolution given 
in equation (36) differs from these earlier, largely empirical, energy crite- 
ria by including the influence of internal work. Albeit this is a lesser 


5ln accordance with equations (14) and (21), the quantities -Bjj/2hh 
denote the rates of change in the thermodynamic displacements (or internal var- 
iables) conjugate to the thermodynamic forces Therefore, the quantity 

- &ijBj^/2hb can be interpreted as the rate of change in the internal work. 
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effect than that of inelastic work, nevertheless, it is believed not to be a 
negligible one, especially when time-dependent effects are present. Addi- 
tional insight is gained from an equivalent expression for the dissipated work 
rate (eq. (34)), that is 


Z 

^if ji 




(G) 


(37) 


which is obtained from equations (29) and (34). The first term in this rela- 
tionship evolves with inelastic deformation and is a measure of fatigue damage. 
(See fig. 2.) The second term evolves with time at internal stress (I 2 is the 
second invariant of internal stress) and is a measure of the interactive creep 
damage. Thus, the interaction between fatigue and creep damage is specified. 
The temperature dependence for this measure of creep damage is accounted for, 
to a large extent, in the thermal recovery function rjj. 


CONCLUDING REMARKS 

A theory of viscoplasticity has been derived from conditions of stability 
and physical arguments, for an initially Isotropic continua that exhibits 
internal damage. This material degradation was incorporated through the 
Kachanov-Rabotnov concept of a net stress. Damage was assumed to be an inter- 
nal variable that evolves isotropically according to a Ponter-Leckie type con- 
stitutive assumption. A potential function was considered that extends the 
Robinson viscoplastic model by including the effects of isotropic hardening 
and material degradation. The yield strength was not considered to be an inde- 
pendent variable; rather, it was assumed to evolve with changes in threshold 
strength and temperature. 

We determined that inelastic strain evolves according to a Prager type 
flow equation, and that Bai ley-Orowan type kinetic equations govern the evolu- 
tion of both internal stress and threshold strength. The internal stress har- 
dens like a Prager hardening rule, whereas the threshold strength hardens with 
dissipated work - not inelastic path length - at a rate that depends on ther- 
mal history. Internal damage was shown to evolve with dissipated work leading 
to a loss of material stability. Dissipated work is the sum of inelastic work 
and internal work. Internal work is a continuum measure of the free energy 
stored in a material due to inelastic deformation. 
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APPENDIX 


This appendix provides the derivations for the flow and evolutionary equa- 
tions given In equations, (28), (29), (33), and (36). 


By using equations, (15), (17), (19), (24), and (25), the flow law 
(eq. (20)) can be written as 


by the chain rule, where 


an 8F 

332 

3E 

vu 

9Snm 

8F aj2 

3£ 

uv 

mn 

3o,j 

an 

K^f(F) 


8T“ 

2n 



3F 

1 



332 

“k2 




(Al) 


(A2) 

(A3) 


8J, 


32 


uv 


uv 


(A4) 


and 


as 


8E 


vu 


3S 


mn 


» S S 
vm un 


nm 


8a, 


U 


= S . 8 . - 8 . . 

ni mj 3 nm ij 


(A5) 


(A6) 


Combining these equations results In 



f(F) E, 


13 


(A7) 


which Is the flow equation (28). 

By utilizing equations (16) to (19) and (24) to (26), the evolutionary 
law for Internal stress (eq. (21)) can be written as 


B 


ij 


-2hjj(G) 


32_3F 

9^2 * 3G 312 


by the chain rule, where 


(A8) 


3E 


vu 



-8 8 
vm un 


(A9) 


8Q 


(AlO) 
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(All) 


and 



dL 


dB. 


mn 


mn 


(A12) 


o|5^j ni mj 3 nm ij 


(A13) 


Combining equations <A2) to (A4) with equations (A8) to (A13), we obtain 

^ij ’ hjj<G)f(F) ^ - 2hjj(G)g(G)B^j (A14) 

which when joined with the flow equation (eq. (A7)) results in 


^ij “ 

where r[)(G) is defined to be 2hj)(G)g<G). This is the evolutionary equation 
for Internal stress given in equation (29). 


We define the following expression: 


W = - K 


^9F 

3F 3K ^ 3G 3KJ 


(A16) 


From equations (18), (19), (25), and (26), we obtain 

3F. .fuiii 

3K- ^3 

2S . . 

3K- ^3 


(A17) 

(A18) 


which when substituted into equation (A16), along with equations (A2) and 
(AlO), gives 


f(F)E,.i: 


H. 


(A19) 


By substituting equations (A7) and (A15) into this relationship, it becomes 




(A20) 


131 


which can be expressed as 







(A21) 


because of equations (15) and (16), and the fact that and are 

deviatoric. This is the rate of dissipated work given in equation (34). 

From equations (24) to (27), the evolutionary law for threshold strength 
(eq. (22)) can be written as 


Z = -h^(Z) 


by the chain rule, where 


and where (from eq. (30)), 

If-. r-J*|fTdt (A24) 

Joining equations (A16) and (A22) to (A24) results in 

^ (r - j‘ i-T «) I - 

where r 2 (Z) is defined to be h 2 (Z)z(Z). This is the evolutionary equation 
for threshold strength given in equation (33). 

By using equations (24) to (27), the evolutionary law for damage 
(eq. (23)) can be written as 


8Q.3G1 ^ 
3F 8K 3G 8K 3Z 


2(Z) 


(A22) 


3Z 


K ^ 
K 3Z 


(A23) 


(0 = D(co) 


3F 3K 3G 3KJ 3u 


by the chain rule, where 


M _ _ K 

3ti) ~ 1 - u 

Combining equations (A16), (A26), and (A27) results in 


(A26) 


(A27) 


CO » D(co) — 

1 - (0 

which is the evolutionary equation for damage given in equation (36). 
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Figure 1 . - In the stress free state fa/h<fa> provides maximum Internal dis- 
sipation during state recovery. 





^ = I (r - ^) d/ 

Figure 2. - Dissipated work In cyclic shear In the absence of thermal recovery. 
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